Sink flow deforms the interface between a viscous liquid and air 

into a tip singularity 
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In our experiment, an interface between a viscous liquid and air is deformed by a sink flow of 
constant flow rate to form a sharp tip. Using a microscope, the interface shape is recorded down 
to a tip size of 1 /im. The curvature at the tip is controlled by the distance h between the tip and 
the sink. As a critical distance h* is approached, the curvature diverges like l/(h — h*) 3 and the 
tip becomes cone-shaped. As the distance to the sink is decreased further, the opening angle of the 
cone vanishes like h 2 . No evidence for air entrainment was found, except when the tip was inside 
the orifice. 

PACS numbers: Valid PACS appear here 



Hydrodynamics very often leads to the spontaneous 
formation of very small structures, such as in the sep- 
aration of a liquid drop from a nozzle 0, , dropcoa- 
lescence EJ, shock waves jj , free-surface cusps j^, , or 
tips 0, H El • Mathematically, all the above phenomena 
correspond to a singularity [T(J, UM of an underlying hy- 
drodynamic equation, in which the hydrodynamic fields 
become non-smooth. Hydrodynamic singularities can be 
classified into two groups. In the first group of 'true' sin- 
gularities the production of small scales is cut off only by 
a microscopic length, such as the size of a molecule. In 
the second group, there exists some hydrodynamic mech- 
anism which is able to 'regularize' the singularity. In the 
past, all singularities which persist for a finite amount 
of time were found to belong to the second group, typ- 
ical regularizing mechanisms being diffusion j^, surface 
tension 0, or the presence of an outer liquid 6]. By 
contrast, in this letter we describe a stationary tip whose 
typical size is not limited by any hydrodynamic mecha- 
nism, such as surface tension. 

Tips were first investigated by G.I. Taylor 0, by plac- 
ing a liquid drop of small viscosity in a viscous flow pro- 
duced by counter-rotating rollers, stretching the drop. 
At its ends, the drop develops tips which become quite 
sharp, yet it has never been established whether the tip 
remains rounded on some scal e Il3j . Tips have recently 
received intense scrutiny 0, Q, 0] owing to their poten- 
tial in producing extremely thin jets , leading to numerous 
microfluidic applications [lJ,[lSll3- I n the viscous with- 
drawal experiment 0, ^| a liquid is withdrawn through 
a straw placed above a liquid-liquid interface, produc- 
ing a sharp "hump". While the curvature of the hump 
appears to diverge as the flow rate is increased, the (ap- 
parent) singularity is actually cut off at a typical size of 
0.2mm 0, and the hump transforms into a jet. 

In the present experiment we replace one of the liq- 
uids by a gas, finding no sign of the disappearance of 
the tip. Hence we are able to track the approach to the 
singularity (i.e. a divergence of the tip curvature) over 
four orders of magnitude, the curvature being described 
by a universal scaling law. Three-dimensional axisym- 




FIG. 1: Sketch of the experimental setup. 



metric flow thus appears to be fundamentally different 
from a two-dimensional one in that surface tension 
does not succeed in regularizing the tip singularity. The 
presence of a singularity also raises exciting prospects to 
produce structures by hydrodynamic forcing whose size 
is only limited by microscopic features [lq . 

Our experimental setup, sketched in Fig.QI consists of 
a27x27xl30 mm cell partially filled with a viscous liquid 
(we used castor oil and two different silicone oils). The 
liquid flows out of a hole (a = 1 and 3 mm in diameter) 
drilled into a circular plate of 0.7 mm thickness. The 
flow rate Q = Q ou t is adjusted to a constant value by 
the ball- valve (1), and measured by weighing. To ensure 
stationary conditions, the flow rate Qi n into the cell is 
adjusted by the valve (2) to be almost equal but smaller 
than Qout- In the course of a typical experiment, the 
cell empties over a period of 2-6 hours at constant Q, 
thus adiabatically changing the filling height H of the 
container. 

All results reported here are for stationary and low 
Reynolds number flow. It thus differs from which 




FIG. 2: Photographs of the air- liquid interface (silicone oil, n 
= 60 Pa.s and Q= 3.9 1(T 3 mLs -1 ). a): h > ft*, b): h = ft*, 
c): h < h* . d): a Gaussian fit to the tip just before the 
singularity, e): the tip as it is about the enter the orifice. 



FIG. 3: Tip curvature C as a function of h for silicone oil 
(77 30 Pa s, 7= 2.13 • 1.0 x 10~ 2 Nm _1 ), a 1 mm, and 
Q=1.7 10" 3 ml s _1 (o), 5.1 • 10" 3 mls- 1 (ffl), 10" 2 mls _1 (•) 
and for a = 3 mm and Q = 3.9 • 10 _2 mls _1 (a). Errors in 
the curvature represent typical deviations between Gaussian 
and parabolic fits, whereas errors of 1.5 % in ft come from the 
finite range of focus. Inset: log- log plot of C versus ft — ft* for 
Q = lO^mls -1 . 



contains a significant time dependence. The interface 
shape (cf. Fig. EJ) is recorded using a microscope with a 
working distance of 20 mm and a a high resolution CCD 
camera. The overall resolution of the optics is 1 pm. The 
benefit of using gravity to drive the flow is that no un- 
steadiness could be detected down to a micron scale, but 
we are limited in the achievable flow rates. Our attempts 
of using a pump for driving failed owing to unsteadiness 
of the interface caused by mechanical vibrations. 

Back-lighting ensures a sharp gray- level gradient at the 
interface, and profiles R(z) were extracted in numerical 
form using image processing software. Axisymmetry is 
preserved to within our experimental resolution. First 
we focus on the development of a sharp tip, shown in 
Fig. 0a), b). To that end we determine the curvature C at 
the tip, by fitting a Gaussian and a parabolic profile (cf. 
Fig. Eld)) to it. Fig.[3|shows the curvature as a function of 
the distance ft from the orifice to the tip for various flow 
rates Q, using the same silicone oil of viscosity n = 30 
Pa s, surface tension 7= 2.13- 10 _2 Nm _1 , and density 
p = 976kg /m 3 . As ft reaches a critical value ft*, the 
curvature diverges like l/(ft — ft*) 3 , as shown in the inset 
for one data set. The same exponent —3 ±0.2 is observed 
for all our data sets. Albeit for a limited range, the data 
of for H is also consistent with ours. 

We now show that the singularity occurs when the 
tip experiences a critical viscous forcing, relative to the 
smoothing effects of surface tension. Namely, in the limit 
of slow flow a point sink in a solid plate produces a ve- 
locity 3Q/(2ttz 2 ) along the line of symmetry p. 140. 
Thus equating viscous forcing by the unperturbed flow 
field with surface tension leads to 3nQ / (27rft* 2 ) ~ 7 for 
the distance ft* where the singularity is expected to oc- 



cur. This is consistent with the data presented in Fig.^J 
if we assume that the point sink is placed a small dis- 
tance h a = 0.29a inside the hole of diameter a = 1mm. 
Thus the singularity is characterized by the local capillary 
number 



Ca(ft) 



reaching a critical value 



3Qr] 



2^(ft±ft a ) 2 7 



Ca(ft*) = Co" = 1.2 ±0.06. 



(1) 



(2) 



This is confirmed in the inset of Fig. fusing all our data 
sets, obtained with two different tube openings a = 1mm 
and a = 3mm. For the shift we used h a — 0.29a, a form 
consistent with a transition from a sink to a Poiseuille 

flow m . 

Utilizing the experimental observation that the pref ac- 
tor of the curvature scales like yfQ for a given liquid, we 
can collapse all the curvature data into a single scaling 
law: 



C= (1.08 ±0.06) 



1/2 



[Qv/(pg)] 
(ft - ft*) 3 



= ii 



lc(ft* + hq) 
(ft - ft*) 3 



(3) 



where £ c = \Ay/ (pg) is the capillary length. Figure [HI 
shows the master curve corresponding to © for all our 
experiments, covering four orders of magnitude in the 
curvature. As expected, © says that the curvature in- 
creases with flow rate, while gravity tends to flatten the 
interface. Surface tension is contained implicitly through 
the dependence on ft*, given by 0. It is instructive to 
also split the r.h.s. of into dependencies on the "in- 
ternal" variable ft and on the "external" capillary length 
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FIG. 4: Critical h — h* with a — 1 mm for silicone oil (77=30 
Pa s) and Q/mls-^1.7 -l(T 3 (o), 5.1 • l(T 3 (ffl), 6.8 • 10~ 3 
(H), 9 • 10" 3 (x), 10" 2 (•), silicone oil ( 77=60 Pa s) and 
Q/mls- 1 ^ 1.4 -10 
and Q/mls _1 = 0.23 
measurement of the orifice position, parameters for castor oil: 
7 = 3.53±0.M0~ 2 Nm"\r7 = 0.93±0.01Pas. Inset: capillary 
number Ca* for all experiments (same symbols, A: a = 3 mm 
with silicone oil (77= 30 Pa s) and Q/mls -1 = 3.9 • 10 -2 .) 
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FIG. 5: Master curve, collapsing all data according to 
Symbols are the same than in Fig. 0] Castor oil data with 
C _1 < 25/ira are not shown, owing to dust perturbing the 
tip. 



£ c , which is required on dimensional grounds. Silicone oil 
and and castor oil have capillary lengths £ c of 1.49 mm 
and 1.94 mm respectively, so we have good experimental 
evidence that it is indeed t c that sets the external length 
scale in rather than the size of the cell. 

We now turn to the interface shape for h < h* , for 
which the tip is singular, and the profile meets the tip 
with a finite slope R f (h). As the hole is approached, 
the slope becomes small, suggesting the use of slender- 
body theory J^I^I^. In this limit, the perturbation 
of the axial flow by the interface is small, so it can be 




FIG. 6: Slope R'(h) versus Ca for all experiments but the 
one with castor oil. Same symbols as in Fig. 0] Line: best fit 
with R'(h) = (0.43 ± 0.1) /Ca. In the shaded region the slope 
is defined as the minimum value of R'(z). 



represented as a line distribution of 2D sources 23]. For 
the slope of the interface at the tip we thus find 



R'(h) = l/[2Ca(h)], 



(4) 



where Ca is given by QJ. Qualitatively, (HJ) results from 
the ratio of the unperturbed flow speed in the z-direction 
and the capillary speed 7/77 contracting the profile in the 
r-direction, while gravity is irrelevant near the tip. 

The slope R' was determined experimentally by fitting 
R(z) with a second order polynomial in \z — h\ and ex- 
trapolating to the tip. The result, shown in Fig.El is in 
excellent agreement with the predicted slope, in partic- 
ular for small opening angles, as expected. The prefac- 
tor is 0.43 ±0.1, only slightly smaller than the predicted 
value of 0.5. To estimate the local capillary number, 
we included the correction h a deduced from h*. This 
we believe to be the first quantitative test of this classi- 
cal slender body theory, previous comparisons being only 
qualitative [24j . The shaded region of Fig. corresponds 
to smooth profiles below the critical capillary number, 
where we took the minimum value of R'(z) as the slope. 
Its value merges smoothly into the singular region. At 
the transition Ca = Ca*, the interface slope assumed a 
universal value of R'(h*) = 0.47 ± 0.06 for all our exper- 
iments. 

Discussion: Our study was motivated by the afore- 
mentioned viscous withdrawal experiment 0, , which 
found the tip to disappear at a typical size of 200 /jul, 
independent of the viscosity ratio A between the two liq- 
uids. This differs markedly from both theory 23] and 
experiment |24| with drops of small viscosity in strong 
flows, which develop tips which are locally similar, but 
which remain stable well into the singular regime. We 
thus expected the disappearance of the tip 8] to be re- 
lated to the far-field boundary condition, which is a flat 
interface for viscous withdrawal, but a curved interface 
for drops. 
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However, at the resolution of the present experiment 
1/im), no air entrainment was ever observed. In fact, 
air bubbles suspended above the sink hole [2^ give re- 
sults very similar to the ones described here, except that 
the curvature scaling does not extend to as small values, 
owing to the bubble curvature. We have considered a 
number of possible causes for the liquid-liquid and the 
liquid-air system to behave differently, yet were not able 
to single out a likely mechanism: 

a) In 15] it was proposed that the minimum tip size 
behaves approximately as C~J n ~ O.l^Ay/ (Kpgj, owing 
to the stabilizing effect of the density difference Ap be- 
tween the two liquids. However, this gives a minimum tip 
size of 150/im for a liquid- air experiment, in disagreement 
with our observations, b) Our viscosity ratios (between 
A = riair I UHquid = 3 • 10 -7 and 2 • 10 -5 ) are smaller than 
those of However, in neither experiment could a sig- 
nificant dependence on A be detected, c) Surfactants are 
known to lead to the ejection of a jet out of the tip of 
bubbles [2(| ("tip streaming"). However, care was taken 
in to differentiate this transient phenomenon from 
the reported instability, d) Finally, it remains to men- 
tion the possible effects of compressibility, solubility, and 
of evaporation into the gas, which distinguish a gas most 
markedly from a liquid. In fact solubility has, in the 
two-dimensional analogue of the present problem, been 
suggested |2g as a mechanism to relieve the pressure that 
builds up inside a cusp, and which leads to its eventual 
disappearance 5]. However, one will require a better 
knowledge of the flow near the tip to estimate the cor- 
rect pressure balance in the present axisymmetric case. 
It should also be noted that once a tip size of 1 micron 
is reached, the mean free path of the air will come into 
play. 

As long as the tip is above the orifice, no air entrain- 



ment occurs. However the indirect evidence of air bub- 
bles appearing at the other end of the tube shows that 
a spout is formed once the tip enters the orifice. In fact, 
this setup has been used to create monodisperse streams 
of bubbles [l(|. When subsequently the filling height 
H is increased again to raise the tip above the orifice, 
no spout is observed. Rather, a stable tip reappears at 
h = 0, giving no evidence for hysteresis. Thus the ex- 
perimental protocol suggested in fails to produce a 
thin spout, a fact that could be related to our far-field 
conditions being different from those of jl8j | . 

A theoretical understanding of the scaling law © re- 
mains a challenge. Likely, such an understanding will 
require the description of the entire interface profile. 
Rescaling the profiles using C _1 as a characteristic scale 
has lead to data collapse only in a very limited region 
around the tip. Evidently the profiles contain several 
length scales that need to be included into a proper de- 
scription. 

In conclusion, we have developed a scaling law for the 
formation of an axisymmetric tip singularity on a liquid- 
air interface. Singular tips are well described by slender- 
body theory. Our analysis raises the possibility of using 
hydrodynamics to produce structures of virtually unlim- 
ited smallness. 
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